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Abstract. We study flocking in one dimension, introducing a lattice model in which 
particles can move either left or right. We find that the model exhibits a continuous 
nonequilibrium phase transition from a condensed phase, in which a single 'flock' 
contains a finite fraction of the particles, to a homogeneous phase; we study the 
transition using numerical finite-size scaling. Surprisingly, in the condensed phase 
the steady state is alternating, with the mean direction of motion of particles reversing 
stochastically on a timescale proportional to the logarithm of the system size. We 
present a simple argument to explain this logarithmic dependence. We argue that the 
reversals are essential to the survival of the condensate. Thus, the discrete directional 
symmetry is not spontaneously broken. 

PACS numbers: 05.70.Fh, 64.60.Cn, 64.60.Ht 

Flocking — the collective motion of a large number of self-propelled entities — is a 
behaviour exhibited by many living beings such as birds, fish and bacteria. Generically, 
flocks are driven, non-equilibrium systems with many degrees of freedom and as such 
they have recently attracted much attention in the physics community [|l], 0, 0, |], |5|, || |7]] . 

Vicsek et al. [|TJ introduced a simple microscopic model in dimensions d > 2 in 
which particles move with a constant speed; they interact only by tending to align with 
their neighbours. Simulations [0] and a continuum theory || showed the existence of a 
low-noise ordered phase in which the mean velocity of the particles is non-zero, i.e., a 
phase in which the rotational symmetry of the model is spontaneously broken. 

Flocking in one dimension (Id) is less relevant to biological systems than in higher 
dimensions, but is nevertheless interesting from a fundamental viewpoint. The models 
studied to date in d > 2 possess a continuous (rotational) symmetry. The Id case is 
necessarily different — since particles are constrained to move either left or right on a 
line, the underlying directional symmetry is discrete. Czirok et al. H have introduced 
an off-lattice model of Id flocking in which particles with continuous velocity variables 
move on a line. Particles tend to move in the same direction as their neighbours but 
'errors' are made due to the presence of noise. Simulations and a continuum theory 
indicate that a continuous phase transition occurs as the noise or particle density is 
varied ||, the ordered phase being characterised by the presence of a large 'flock'. 
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When a phase transition occurs from a disordered to an ordered phase, it 
is usually accompanied by spontaneous symmetry breaking; in the ordered phase, 
the spontaneously broken symmetry is accompanied by ergodicity breaking in the 
thermodynamic limit. In an ordered system with discrete symmetry, such as the Ising 
model in d > 2, one expects the 'flip time' — the time the system takes to move from 
one symmetry-related state to another — to diverge exponentially with system size. 

In this work we study a simple lattice model of flocking in Id. We show, using finite- 
size scaling, that the model has a continuous phase transition from a homogeneous 
to a condensed phase. However, the condensed phase is not symmetry-broken, but 
alternating — we find that the mean particle velocity alternates its sign on a timescale 
which grows only logarithmically with system size; moreover, these reversals are essential 
to the maintenance of the order. 

We now define the model that we study. We consider N particles on a periodic Id 
lattice of L sites. The particle density is p = N/L. Each particle, labelled by /i, has a 
position £ {1,2,..., L} and a velocity (or direction) t> M = ±1. To update the system, 
a particle is chosen at random. Its velocity t> M is flipped with probability and then 
Xfj, — * %[i ~\~ Vp. The flip probability is given by 

W, = [1 - (1 - 2rj)v fl U(x^} /2, (1) 

where U(y) is the velocity of the majority of the particles (including particle /x) at site 
y and its two nearest neighbours; we take U = when there is no majority. In brief, 
particles acquire the velocity of the majority of their neighbours with probability 1 — rj. 

While our lattice model is somewhat simpler than the off-lattice model of reference 
f|, it is conceptually similar. However, whereas the off-lattice model has only one 
source of noise (a random perturbation added to the velocity of a particle at each time 
step), the noise in our model has two distinct sources. The first is due to the flipping 
of particle velocities and its strength is parameterised by 77; the second is due to the 
random sequential dynamics of the model. Thus, unlike the off-lattice model, our model 
does not become deterministic in the limit 77 — > 0. 

Before showing that the model exhibits a continuous phase transition, we first 
discuss the qualitative behaviour. When 77 is large, the system has a fairly homogeneous 
steady state, illustrated in figure |l|(a) by a space-time plot for 77 = 0.2. However, when 77 
is small, 'condensation' occurs — a large fraction of the particles are contained in a single 
large flock, which we refer to as the condensate. A space-time plot of the steady state 
for 77 = 0.02 is shown in figure [j](b). The direction of motion of the condensate reverses 
at fairly regular intervals. These reversals are sometimes, but not always, triggered 
by a collision with a smaller flock travelling in the opposite direction. Just after the 
condensate has reversed its direction, it is very dense. It then gradually spreads out, 
becoming more diffuse at the front than at the back, before reversing again. The steady 
state clearly has no time-reversal symmetry and is alternating, by which we mean that 
the system reverses its mean velocity at stochastic intervals. A quantity measuring 
the spatial order in the system, such as the variance of the number of particles per site, 
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Figure 1. Space-time plots of steady-state systems with (a) 77 = 0.2 and (b) 77 — 0.02. 
In both cases, p = 1 and L = 1000. A darker grey level indicates a higher particle 
density. There is one timestep between each snapshot in (a) and 5 in (b). 




Figure 2. (a) Log-linear plot of (T) and (T 2 ) 1 / 2 as a function of L for 77 = 0.02 and 
p = 1. The straight lines are fits to the data, (b) Log- log plot of (T) and (T 2 ) 1 / 2 as a 
function of 77 for L = 4000 and p = 1. The solid lines are fits to the data for 77 < 0.02 
with exponents of —0.358 ± 0.01 and —0.366 ± 0.01 respectively. 



decreases gradually as the condensate spreads, then increases suddenly as the condensate 
'flips', before decaying again. 

If a condensate were never to reverse its direction of motion, it would eventually 
become so diffuse that it would cease to exist — the alternating character of the steady 
state is essential for the existence of the condensate. The condensate has no preferred 
direction of motion and we conclude that directional symmetry is not broken. 

We now discuss the behaviour of the condensate in more detail. Let p{T) be the 
probability that the condensate reverses its direction (i.e. that the mean particle velocity 
changes sign) a time T after the previous reversal. Figure 0(a) shows a log-linear plot of 
(T) and (T 2 ) 1 / 2 as a function of system size L for fixed 77 = 0.02. Both are proportional 
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Figure 3. (a) Log-log plot of w(r)) for p = 1. The solid line is a fit to the data, (b) 
Linear-log plot of the scaled reversal time distribution (T)p(T) as a function of scaled 
time T I (T) for p = 1 and various values of L and 77. 



to logL. Thus, while no spontaneous symmetry breaking occurs, the mean time between 
reversals does diverge weakly in the thermodynamic limit, as does the variance. 

Figure @(b) shows a log-log plot of (T) and (T 2 ) 1//2 as a function of 77 for fixed 
L = 4000. For 77 less than about 0.02, both (T) and (T 2 ) 1 / 2 diverge as power laws in 77. 
The data suggests that the exponent, which we call A, is the same in both cases with a 
value of —0.36 ± 0.01. We have seen in figure 0(a) that (T) oc logL for fixed 77 and L 
large; we denote the constant of proportionality by -1/7(77). Figure |3](a) shows a log- log 
plot of w{rj) and a least-squares fit to the data gives an estimate of —0.34 ± 0.02 for the 
slope; it seems likely that this is equal to A and so we conclude that (T) ~ r] x \ogL for 
large L and small 77. 

The fact that (T) and (T 2 ) 1 / 2 are both proportional to logL for fixed 77, and 
both proportional to r] x for fixed L suggests that the distribution of reversal times 
may have the scaling form p(T) ~ f(T/(T))/(T) asymptotically for large L and small 
77. Figure |3](b) shows the scaled distribution for various values of L and 77 — the data 
collapse is fairly good. 

We now present a simple argument which may explain why a condensate containing 
N particles flips direction in a time of order log N (or equivalently log L if the density is 
fixed) for large N. Consider a system in which, at t — 0, all N particles are positioned 
at x = and have velocity +1. Let q(x, t) be the probability that a given particle in this 
initial flock has position x at time t. In the limit 77 — > 0, each particle hops forward with 
probability 1/N in a time At so that in the limit iV — > 00, At = 1/N, the distribution 
q(x, t) is Poisson: 



For 77 nonzero, this distribution will be modified by the flipping of particle velocities. 
For example, particles will leave the back of the condensate when 77 is nonzero but here 
we neglect this effect. 
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We now ask the following question: if some perturbation or fluctuation (due to 
collision with a small flock or the spontaneous flipping of particles at the front of the 
condensate) should cause all the particles in the condensate having position x > z to 
flip velocity, will this in turn cause the particles with position x = z to flip? We make 
the approximation that the particles with x > z form a 'shock' moving back through 
the condensate so that at time t they all have position z + 1 and velocity —1. Then the 
particles with position z will flip if R(z,t) > 1, where 

1 00 t 00 z! t x 

R(z,t) = - — V q{x,t) = — — V -. rr . (3) 

K,J q(z-l,t) + q(z,t) x ^ +i qK ' z + t^[(z + x)\ K) 

For a given position z we define a critical time t c (z) through the solution of R(z, t c ) = 1 
so that for t > t c (z), the particles at position z are susceptible to flipping. 

It is easy to see that R(z, t + 1) > R(z + 1, t) so that if the particles at z + 1 are 
susceptible to flipping (i.e. if R(z + l,t) > 1) then those at z are subsequently also 
susceptible. Hence, the reversal time of the entire condensate is determined by the time 
one must wait for the leading particles to flip and it suffices to find t c (z) in the limit 
z — ► oo. It is straightforward to show that R(z, t) — ► s 2 / (1 — s 2 ) as z — > oo with t = sz; 
this implies that t c (z) — > z/y/2. Therefore, £ c is approximately proportional to z for 
large z. 

For a condensate containing a large but finite number of particles N, we estimate 
z*(N,t), the position of the leading edge of the condensate, through the expression 
q(z*,t) = 1/N. For the Poisson distribution when t oc z* we find that z* ~ logiV. 
Therefore our estimate for the time at which the condensate becomes susceptible to 
reversal is t c (z*) which, being proportional to z*, is of the order of log N. 

Since this simple analysis makes the assumption that one can decouple the spreading 
of the condensate and the fluctuations which may cause it to reverse its direction, it 
makes no predictions about the non-trivial r] dependence of the reversal time and the 
scaling behaviour of the distribution of reversal times. 

This completes our discussion of the condensed phase; we now turn to a numerical 
finite-size scaling analysis || of the phase transition between the condensed phase and 
the homogeneous phase for fixed p — 1. We have performed Monte Carlo simulations 
for system sizes between L = 50 and L = 2000, averaging over between 5 x 10 6 and 
2 x 10 7 time steps in the steady state for each set of parameters. We have found that 
the absolute value of the mean particle velocity V, defined by 

(4) 

is a convenient order parameter since, between reversals in the condensed phase, the 
majority of particles have the same velocity, that of the condensate. During the long 
(but O(logL)) time intervals between reversals, V fluctuates about some well-defined 
mean value. Figure f|(a) shows a plot of (V) against r\ for several different system 
sizes. The crossover from the disordered to the ordered regime becomes sharper with 
increasing system size, suggesting the presence of a continuous phase transition. 



N 
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Figure 4. (a) (V) as a function of rj for various system sizes with p = 1. (b) Ul(t]) 
for various system sizes in the vicinity of rj c , given by the common intersection point. 
The solid lines are cubic polynomial fits to the data. 



In calculating critical quantities, it is useful to define both the 'static susceptibility' 
and the fourth-order cumulant given respectively by 

X (L) = L[(V 2 ) L - [V)l] ; U L = 1- (5) 

with the angle brackets indicating steady state time averages. In analogy with magnetic 
systems we assume that, near criticality, (V) ~ ij] c — rj)^, x ~ \Vc~ ??|~ 7 and £ ~ \r] c — ri\~ u , 
where £ is the correlation length. 

Standard finite-size scaling theory || leads to the following relations at criticality: 

(V) L ~L-^; X (L)~L^; U L = U*. (6) 

The critical point rj c can therefore be identified as the value of 77 at which Ul takes on 
its fixed point value U*. Figure f|(b) shows that the curves Ul{t]) do indeed have a 
common intersection point at 7] c = 0.1325 ± 0.001. Also, we find U* = 0.30 ± 0.01. 

The slopes of the functions Ui^q) at rj = i] c can be used to calculate the critical 
exponent v since || 

l/V = log (dUbL/dUL)\ Vc /log b + corrections to scaling. (7) 

We have found that Ul(j]) is approximately linear near r\ c for L less than about 1000, 
allowing a reasonably precise estimate of the derivative in (0). In figure |5|(a) we plot 
estimates of v obtained using (^) against l/log& for various L. Extrapolating b — ► 00, 
thus taking account of corrections to scaling, we estimate v = 2.57 ± 0.05. 
The other independent exponent 7/V is given by 

7/V = \og[x(bL, i] c )/x(L, ^ c )]/log6 + corrections to scaling. (8) 

Figure |](b) shows the resulting estimates of 7/V plotted against l/log& for different 
values of L. As before, we extrapolate b — ► 00 and estimate 7/V = 0.451 ± 0.005. 
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Figure 5. Estimates of (a) v and (b) 7/1/, obtained using (Q) and (|J), plotted against 
l/log& for various system sizes. The lines in both are guides to the eye. 

The finite-size scaling formalism relies on the fact that the hyperscaling relation 
d — y/v + 2(3/v holds [|J. At r) = r\ c one can calculate (3 /v via 

(3/v = — \og[(V) bL / (V) l]/ logb + corrections to scaling. (9) 

We were unable to estimate (3/v as precisely as y/v but, within error bounds, we have 
found that hyperscaling does indeed hold. Thus, our final estimates for the critical 
exponents are (3 = 0.705 ± 0.02, 7 = 1.16 ± 0.04 and v = 2.57 ± 0.05. 

In reference [[| (3 was estimated (without finite-size scaling) for an off-lattice 
model of Id flocking to be 0.6 ± 0.05. While our value of /3 = 0.705 ± 0.02 is not 
in close agreement with this, it is still possible that both models are in the same 
universality class. Measurement of U* in the off-lattice model would provide a good 
test of universality. 

Above, we have used standard finite-size scaling to study the order parameter 
distribution. The success of this approach is surprising when one considers that the 
order parameter V does not fully reflect the alternating nature of the steady state in 
the condensed phase — the 'spreading' of the condensate is not captured. 

We have argued that spontaneous symmetry breaking does not occur in our model 
and that the reversal times are O(logL). This contrasts with equilibrium systems where 
reversal times exponentially large in L result from surmounting free energy barriers. It 
is difficult to see how O(logL) reversal times could be explained in an analogous fashion; 
the reversal mechanism is fundamentally nonequilibrium. 

Finally, we comment on the effects of imposing a restriction on the local particle 
density in our model; this is achieved by setting the maximum number of particles 
allowed on a site to be M. Not surprisingly, we have found that for any finite M < N, 
the condensed phase is suppressed; instead, for small rj, the system forms 'domains' of 
maximally occupied sites. Each of these domains comprises two competing sub-domains 
composed of particles with velocities +1 and —1. The interface between the two sub- 
domains performs a random walk as a result of the flipping of particle velocities and the 
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dynamics resembles that of a Id ferromagnet. 



Acknowledgments 

We thank M E Cates for many useful discussions and a critical reading of the manuscript 
and N B Wilding for helpful advice on finite-size scaling. 



References 



Vicsek T, Czirok A, Ben- Jacob E, Cohen I and Shochet O 1995 Phys. Rev. Lett. 75 1226 
Czirok A, Stanley H E and Vicsek T 1997 J. Phys. A: Math. Gen. 30 1375 
Toner J and Tu Y 1998 Phys. Rev. E 58 4828 

Bussemaker H J, Deutsch A and Geigant E 1997 Phys. Rev. Lett. 78 5018 



Czirok A, Barabasi A and Vicsek T 1997 Preprint fcond-mat/9712154 
Albano E V 1996 Phys. Rev. Lett. 77 2129 

Shimoyama N, Sugawara K, Mizuguchi T, Hayakawa Y and Sano M 1996 Phys. Rev. Lett. 76 3870 
Binder K and Heermann D W 1997 Monte Carlo Simulation in Statistical Physics. Solid-State 
Sciences (Berlin Heidelberg: Springer), 3rd edn. 



